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It is shown that statistical mechanics is applicable to quantum systems with finite 
numbers of particles, such as complex atoms, atomic clusters, etc., where the resid- 
ual two-body interaction is sufficiently strong. This interaction mixes the unperturbed 
shell-model basis states and produces "chaotic" many-body eigenstates. As a result, an 
interaction-induced equilibrium emerges in the system, and temperature can be intro- 
duced. However, the interaction between the particles and their finite number can lead 
to prominent deviations of the equilibrium occupation numbers distribution from the 
Fermi-Dirac shape. For example, this takes place in the cerium atom with four valence 
electrons, which was used to compare the theory with realistic numerical calculations. 

1. Introduction. Statistical behaviour is usually established in the limit of a 
large number of particles n. Moreover, quantitative results can be obtained if cor- 
relations between the particles are somehow weak. This means that the interaction 
between the particles is to be neglected, or, a more realistic possibility, an appropri- 
ate mean field theory is chosen. The latter results in the picture of free quasiparticles 
moving in the effective self-consistent field created by the constituents. 

In this limit the temperature T is a well-defined physical quantity and all equi- 
librium characteristics can be found. For example, for a gas of noninteracting 
fermions this results in the famous Fermi-Dirac distribution (FDD) of the occupa- 
tion numbers: n a = l/(exp[(e a — ji)/T] + 1) . There are many real complex systems, 
such as compound nuclei, rare-earth atoms, molecules, quantum dots, etc., which 
do not satisfy the conditions for FDD to hold. The number of active particles in 
these systems can be relatively small (~ 10), and the interaction between them 
(even the residual interaction in the mean-field basis) is large, i.e., greater than the 
energy intervals between unperturbed basis states. However, this interaction makes 
up for the absenccj of a heat bath, and promotes the onset of "randomization" and 
quantum chaosiratl, which gives one a possibility to talk about some kind of equilib- 
rium in the system, and pursue the development of a statistical theory for few-body 
Fermi systems 0. 

2. The Ce atom. The cerium atom has one of the most complicated spectra 
in the periodic table. Its electronic structure consists of a Xe-like Is 2 . . . 5p 6 core 
and four valence electrons. A large difference in the energy scales of the core and 
valence electrons allows us to neglect excitations from the core and consider the 
wave function of the core as a "vacuum" state |0). Accordingly, the four active 
electrons form the spectrum of Ce below the ionization threshold. 

The calculations are performed using the Hartree-Fock-Dirac (HFD) and con- 
figuration interaction (CI) methods (seetlfor details). A self-consistent HFD calcu- 
lation determines the mean-field potential, which is then used to calculate the basis 
set of single-particle ortho-normalized relativistic states |a) = \nljj z ) with energies 
e a . This procedure defines the zeroth-order Hamiltonian of the system, ~ 

£ Q a l a ct ■ The unperturbed multi-particle basis states |fc) constructed from the 
single-particle states, |fc) = a^a^a^a^ |0), are eigenstates of H^: H^\k) = 
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E^\k), where E^ — ^) E Q Ha' and = (fc|a^a a |fe) are the occupation num- 
bers equal to or 1. 

The total Hamiltonian H of the active electrons is the sum of the zeroth-order 
mean-field Hamiltonian and the 2-body residual interaction is 
V — \ X)a/37(5 ^a^7i5 a Q a /5 a 7 a i5 ■ This interaction contributes to the diagonal and off- 
diagonal matrix elements between the multi-particle states |fc). The diagonal part 
represents an additional energy shift of the state k, and can still be expressed in 
terms of the occupation numbers. The off-diagonal matrix elements are responsible 
for mixing of the multi-particle basis states. 

Complete diagonalization of the operator H = H4 + V within the space of the 
\k) basis states produces "exact" energies Ei and stationary states 

H\i) = E i \i), |i)=£<7f|fe>. (1) 
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We included 14 relativistic subshells nlj in the calculation (6s, 7s, 6p, 7p, 5d, 
6d, 4/, and 5/), and performed exact diagonalization of the NxN Hamiltonian 
matrix in a Hilbert space with JV ~ 8 x 10 3 , obtained by truncating the complete 
set of the shell-model atomic configurations. To subtract additional symmetries 
our numerical results are obtained for the even states of Ce with the total angular 
momentum projection set to J z = 0. 

The relativistic atomic subshells nlj are g s — 2j + l degenerate, therefore, we 
consider average occupation numbers for the subshells s: h s = gs 1 J2 a es a a a a- 
When the number of active particles is small, the occupation numbers for any 
eigenstate ris = (i\n s \i) show strong level-to-level fluctuations, and it is more 
instructive to look at spectrally averaged values 
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n e (E) = {i\n s \i) = > \Cl l> (k\h s \k) , (2) 



where the overline means averaging over the eigenstates i within some energy in- 
terval around energy E. 

A typical distribution of the occupation numbers calculated at the excitation 
energy of 3.75 eV above the atomic ground state is shown in Fig. [l] as a function 
of the single-particle energy e s of the orbitalal. One can see that the distribution 
does not look at all like a monotonically decreasing FDD, as some higher-energy 
orbitals have larger n s (E) than the lower ones. Moreover, this behaviour persists 
over the whole energy interval from the ground state to the ionization potential. 
For example, the lowest even state of Ce has a configuration of 4/ 2 6s 2 , while the 
FDD would tell us that all 4 electrons must be placed in the lowest 4/-orbital, when 
the energy of the system is low. At first sight such a strong deviation from the FDD 
in a strongly interacting Fermi system speaks against any possibility of statistical 
description of the system. However, we show that strongly interacting orbitals can 
be properly incorporated in the canonical ensemble description of the system, and 
thermally averaged occupation numbers n s (T) derived. 

3. Statistical model. Let us perform a statistical calculation of the occupa- 
tion numbers for a system of n particles distributed over r orbitals with energies e s 
and degeneracies g s (s = 1, . . . , r). We will assume that the two-body interaction of 
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Figure 1: Occupation numbers n a {E) (see Eq. (^)), for the even states of Ce at the excitation 
energy of E = 3.75 eV versus the single-particle energies e s of the orbitals. 

any two particles in the orbitals s and p is J7 sp , where both the direct and exchange 
terms are included: 
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sp gs(g P -S sp ) 



The energy of a particular many-particle state fc is now given by 



E k 



E-^+EE 



N S {M P -5 S 

1 + 6 S r, 



-u, 



(3) 



(4) 



where TVs is an integer number of particles in the orbital s (0 < A/" s < g s ), and 
^ s A/" s = ft. The state k is defined by specifying the orbital occupation numbers 

J\f 8 , and is G^-degenerate, where Gk — ]X=i )• 

In the canonical ensemble at temperature T the probability of finding the system 
in the state fc is given by 



w k = Z- x G k exp(-E k /T) , where Z = V G k exp{-E k /T) 
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(5) 



and the sum over k runs over all possible multi-particle states (possibly, with the re- 
striction of parity). The average occupation numbers n s (T) = J\f s /g s are calculated 
as 



n s {T) = g\ 



r(k) 

where Ms is the number of particles in the orbital s in the multi-particle state 
k. The energy of the system at a given temperature is E(T) = J2k EkWk- If we 
know E(T), we can find inverse function T(_E),and use it in (6) in order to compare 
with the data from (||). This comparison is shown in figure 2 and good agreement 
is observed. 

4. Discussion. This agreement means that a kind of equilibrium is indeed in- 
duced in the system due to the interaction between the particles ( "micro-canonical" 
distribution). Moreover, the averaging over it yields results close to those over a 
canonical ensemble (||), with the temperature chosen to reproduce the total energy 
of the system. This equivalence is always true for large systems where any, albeit 
weak, interaction between particles leads to equilibrium. However, in a few-particle 
system the residual two-body interaction must be strong to produce chaotic eigen- 
states and facilitate statistical description. 
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Figure 2: Comparison of the orbital occupancies g s n s (E) obtained from the exact diagonahzation 
(solid and dash-dot lines) with g s n s (T(E)) obtained from our statistical theory (dotted lines). 

Note, that although the temperature-based description is valid for our 4-particle 
system, the orbital occupancies could not be described by the FDD (Fig. [I]). The 
FDD is inapplicable to our system because of the strong interaction between the 
particles [second term on the the right-hand side of Eq. However, the deviation 
from the FDD is determined not by the magnitude of U sp , but rather by the size of 
their fluctuations. To see this assume for a moment that U sp = U are the same. In 
this case the double sum in Eq. ([|) just shifts all energies by ^-N(N — 1), and the 
statistical properties of the system are the same as for noninteracting particles. If 
U sp are different for different orbital pairs sp one can still introduce some average 
interaction U and subtract this "background" interaction from the interaction term 
in Eq. (||). This procedure will effectively suppress the interaction term, since the 
summands in expressions like J2 s < P (Us P — U)N P have different signs. Note that the 
introduction of (energy-dependent) U is equivalent to a mean-field approximation. 
This approximation is good if the fluctuations of U sp from one orbital to another 
are relatively small. This formulates a condition for the FDD to be valid. In the Ce 
atom the situation is just opposite. The 4/ orbital has a much smaller radius than 
any other orbital, hence the Coulomb interactions U^f^f or U^f s are much greater 
than any other U sp (here s and p are orbitals other than 4/). 
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4. Of course, when this interaction is strong one needs to introduce some new 
mean-field orbital energies e s that would incorporate the effect of such inter- 
action. The value e s will inevitably depend on the distribution of the other 
electrons, and hence on the excitation energy of the system. However, even 
when s s are used there is a noticeable deviation from the simple FDD distri- 
bution. 
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